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We predict that a Majorana quantum wire, driven into a trivial phase by tilting the external mag- 
netic field as observed in experiment [V. Mourik et al., Science 336, 1003 (2012)], can be restored 
into the topological phase by applying a supercurrent in the proximate superconductor (SC). A 
tilted external magnetic field breaks spatial inversion symmetry and closes the superconducting gap 
in the wire when the tilt angle exceeds a critical value, leading to a topological phase transition. In- 
terestingly, we show that a supercurrent applied in the SC is able to restore the Majorana end-states. 
Using Abelian bosonization, we further confirm this result in the presence of electron-electron inter- 
actions and show the analogy of this phenomenon to the commensurate-incommensurate transition 
driven by doping a one-dimensional Mott-insulator system. The present results have important 
applications in e.g., realizing the supercurrent assisted braiding (SAB) of Majorana fermions, which 
proves highly useful in topological quantum computation with realistic Majorana networks. 

PACS numbers: 71.10.Pm, 74.45. +c, 03.67.Lx 



The study of topological superconductors (SCs) which 
host Majorana zero bound states (MZBS) has developed 
into a remarkably lively and rapidly growing branch of 
condensed matter physics, driven both by the pursuit of 
exotic fundamental physics and the applications in fault- 
tolerant topological quantum computation (TQC) [IH3]- 
A MZBS is a self-Hermitian and charge-neutral excita- 
tion which is predicted to exist in the vortex core of a two- 
dimensional (2D) (p + ip)-w&ve SC [4], and at the edges 
of a one-dimensional (ID) p-wave SC [5J [BJ. However, in- 
trinsic p-wave superconductivity is not necessary to ob- 
serve MZBS: recent proposals have shown the equivalence 
of topological insulator/ s- wave SC heterostructures [7HS] 
and spin-orbit (SO) coupled semiconductor/ s- wave SC 
heterostructures with Zeeman splitting 10 16J to p-wave 
SCs. In such devices, the SO interaction drives the orig- 
inal s-wave SC into an effective p-wave SC, leading to 
MZBS in the case of odd number of subbands crossing 
the Fermi energy. It has been predicted that MZBS at 
the edges of a ID topological SC are experimentally de- 
tectable in the differential tunneling conductance dl/dV 
at the interface with a normal contact, via the emer- 
gence of a zero-bias peak (ZBP) of height 2e 2 /h (at zero 
temperature) [174420J . The recent experiment using InSb 
nanowires (NWs) connected to normal leads via tunnel 
barriers and in proximity to s-wave SC NbTiN has shown 
a suggestive ZBP in the dl/dV spectra, which disappears 
upon tilting the external magnetic field, in accordance to 
theory 21 . The height of the ZBP is however much less 
than 2e 2 //i, a fact that may be due to thermal broadening 
or multiband effects [12 [52]. ^ more rigorous verification 
of MZBS is at present under intensive investigation. 

Motivated by these recent findings, in this work we 
study a Majorana quantum NW driven into the trivial 
phase by the application of a tilted external magnetic 



field which breaks ID spatial inversion symmetry (SIS) 
[2"U] , as observed in the experiment [5T] . Interestingly, we 
predict that the NW can be restored into the topological 
phase by applying a supercurrent in the s-wave SC. We 
show that a supercurrent with proper magnitude and di- 
rection can compensate for the detrimental effects of the 
tilted magnetic field, therefore restoring the topological 
phase. Our result provides both useful experimental tests 
for topological superconductivity in the lab, and impor- 
tant practical applications in TQC via a supercurrcnt- 
assisted braiding (SAB) scheme for MZBS in realistic 
Majorana networks, where NW segments are generally 
tilted from the direction of the external magnetic field. 

We start from the model of a ID SO-coupled NW in 
proximity to an s-wave SC, with a Zeeman splitting V x 
given by an external magnetic field parallel to the NW. 
A topological phase transition from trivial to topologi- 
cal superconducting phase occurs by tuning V x beyond a 
critical value V c = (p 2 + lA^ 2 ) 1 / 2 pH HO [23], where 
fi and A s are chemical potential and induced s-wave 
SC order parameter in the NW, respectively. Allow- 
ing a tilt angle 6 between the magnetic field and the 
NW direction leads to the angle-dependent Zeeman field 
{V x ,V y ) = Vq (cos#, sin#). The Hamiltonian of the sys- 
tem is therefore given by H = Hq + H s , where 



d 2 
2m* 



dx c l ( x ) 

+ {i^Rd x + V y ) (T y + V x a x \ a a , c a > (x) , (1) 

H s = j dx [A s c T (x) c+ (z) + H.c] , (2) 

with m* = 0.015m e the effective mass of electrons in 
the InSb NW [21 ], X R the Rashba SO coupling coeffi- 
cient, and c a (x) the electron annihilation operator for 
the spin state a. The term V y a y , occurring due to a 
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finite tilt-angle 9, breaks SIS of the NW [5U]. This 
can be seen directly in the Hamiltonian under the ID 
space-inversion transformation x — > —x, (y,z) —¥ (y,z), 
which leads to (k x ,a y ) — > (—k x ,~a y ), and a x — » a x . 
The broken SIS leads to an asymmetric single-electron 
dispersion relation ^ e^' 



for 



for ifn, where eS ± ' 1 = 



k x /2m*±^V x 2 + (V y - X R k x ) . Accordingly, the Bogoli- 
ubov quasiparticle spectra for uniform A s also become 
asymmetric E (k x ) ^ E{—k x ) (Fig. [I] (a)). In particular, 
when 9 is greater than a critical value 9 C (Vq, A s , /i), the 
minimum (maximum) energy of the electron-like (hole- 
like) states becomes negative (positive). Therefore the 
bulk gap of the NW closes (red dashed curves in Fig. [I] 
(a)). This leads to a topological phase transition at 
9 = 9 C , and for 9 > 9 C the system is a trivial SC. 
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FIG. 1: (a) Vanishing of the bulk gap by increasing tilt angle 
9; (b) Restoring of the bulk gap at 6 — 0.2n by applying su- 
percurrents; (c) Phase boundary between the topological su- 
perconducting (T.S.) and trivial phases with a supercurrent; 
(d) Superconducting bulk gap versus optimal supercurrent. 
Other parameters are taken as Vo = l.OmeV, A s = 0.5meV, 
p — 0, and SO coupling energy E BO = m*\ 2 R — O.lmeV for (a- 
d). The critical angle in this parameter regime 9 C w 0.1687T. 

We proceed to show that the topological phase can 
be restored at 9 > 9 C by a supercurrent applied in 
the proximate SC. In the presence of a uniform su- 
percurrent, the induced SC order parameter acquires a 
position-dependent phase A s (x) = |A s |e^ (;r) , related fermions c L/R 



E k [4 } ~lA c l~ck,- + \ Ek (A.e i »c fc) _a_ Jt _ a ,_ +H.c 

where \k = tan -1 [(V^ — \ R k) /V x ] and e 
k > and < 9 < it. In the absence of supercur- 
rent, electron states with opposite momenta ±fe are off- 
resonant and the formation of Cooper-pairs with zero 
center-of-mass momentum is weakened. For a supercur- 
rent applied along +x direction (i.e., a > 0), the pairing 
Hamiltonian pairs up states with momenta k and — k — a 
which are closer in energy, favoring the formation of a 
Cooper pair with center-of-mass momentum a. The ap- 
plication of a supercurrent therefore allows to compen- 
sate for the band asymmetry induced by the tilted mag- 
netic field, strengthening the superconducting bulk gap. 

In Fig. [l] (b) we show that the bulk gap, which for 
J s = vanishes at 9 C = 0.2-7T, reopens with the applica- 
tion of a supercurrent in the +x direction, and attains 
its maximum at the optimal value J s = J° p (green solid 
line). Further increasing J s suppresses again the bulk 
gap (black dotted line) . Physically, in this regime the su- 
percurrent overcompensates the effect of V y and induces 
again an off-resonant situation. Fig. [T] (c) shows numer- 
ical results for the phase diagram as a function of J s 
and 9, with the coherence length in NbTi SCs taking the 
typical value £ < lOnm [25 . The blue curve represents 
the optimal supercurrent J° p (9) , and the red curves give 
the boundaries of the topological and trivial phases. For 
J s = 0, the phase becomes trivial when 9 C < \9\ < tt/2, 
while applying a supercurrent along +x (or —x, depend- 
ing on the sign of 9) can restore the topological phase 
(Fig. gc)). In contrast, for 9 = the NW has SIS, and 
applying a supercurrent may destabilize the topological 
phase [35]. The bulk gap versus optimal current J° p is 
given in Fig. [I] (d) , from which one finds a vanishing bulk 
gap Eg ( J° p ) only for 9 = 7r/2,37r/2, indicating that the 
MZBS can always be restored by the supercurrent unless 
the magnetic field is perpendicular to NW. 

To determine if the above results are robust against 
the presence of electron-electron (e-e) interactions in the 
NW, we introduce here the Abelian bosonization frame- 
work. At low energies, linearization of the dispersion re- 
lation e, around the Fermi energy E R generates asym- 
metric left (right) Fermi momenta k^ (k R ) and left (right) 
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Fermi velocities i>l(_r) = h 



\dks y k 
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as a conse- 



quence of the broken SIS. We next introduce the stan- 
dard bosonic representation of left- and right-moving 
-i(=FV3-i?) w ]j ere t ne bosonic fields 
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to the supercurrent in SC through the relation J s = 
2n s eha[l — (a£) 2 ]/m e [53], with a = V<f>(x) a uniform 
phase-gradient, and m e , n s , and £ the electron mass, 
superconducting carrier density and coherence length, 
respectively. The applied supercurrent density is re- 
quired to be less than the superconducting critical cur- 
rent J c = 4n s eh/ (3v3m e £) [53]. The physics of the 
problem can be seen more transparently by projecting 
the Hamiltonian onto the lower subband H w = 



ip (x) , $ (x) are conjugate canonical variables obeying the 
commutation relation [tp (x) , $ (y)} — «7rsign (y — x) /2 
and a ~ fc^ 1 is the short-distance cutoff of the continuum 
theory [57] . Physically, the field ip (x) represents slowly- 
varying fluctuations in the electronic density bp (x) = 
— d x <p (x) /it, and *& (x) is related to the phase of the SC 
order parameter through the relation — ic (x) d x c (x) cx 
cr (x) cl (x) cx e* 2 *( x ), where we have neglected less rele- 
vant higher-order terms in d x d (x). A short-range inter- 
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action Hi nt — g J dx c R (x) cr (x) c\ (x) cl (x) acquires a 
simple form in bosonic representation and the low energy 
Hamiltonian is given by 



H 



dx 



^ (d x 0) 2 + 

2tt V x ' 2irK 



(d x ipf + —d x ^d x -d 

7T 



+ sin [2tf (x) + (a- 5k F ) x] 
wo, 



(3) 



where interactions are encod ed in the dimc nsionless Lut- 
tinger parameter K = ^/ jfj^g^ j [with K < 
1 (K > 1) for repulsive (attractive) e-e interactions], 
u = (vji + Vl) K/2, and the effective p-wave SC or- 
der A p = Assin^^ 1 — x-k R - a \ ^ rp-^g (jijuengioniess p a _ 

rameter rj — (vr —vl)/ {vr + vl) quantifies the band- 
asymmetry and finally Skp = &l — For A p = 
and rj <C 1, the above model describes an asymmetric 
Luttinger liquid (LL) fixed-point, with broken SIS and 
right- and left-going ID plasmon excitations traveling at 
different velocities vl ^ vr 28 . In the presence of a 
small A p , the low-energy properties of the model are cap- 
tured by the perturbative renormalization-group (PRG) 
approach around the LL fixed-point [2"§l - f3"2"] . For simplic- 
ity, in what follows we neglect the marginal symmetry- 
breaking term ~ rjd x Lpd x 'd assuming that r\ <C 1. Imple- 
menting a standard PRG procedure that leaves invari- 
ant the LL Gaussian fixed-point under the change in the 
short-distance cutoff a (£) = a$e l — > a (£ + d£) allows to 
obtain the RG-flow equations: dK/dl = y 2 Jq (Spa (£)) , 
dy/d£ = (2-K-^y and du/d£ = -y 2 uKJ 2 (5pa(£)), 
with Sp = a — Skp [27J. Here J n (z) is the n-th or- 
der Bessel function of the first kind and y = A p ao/u is 
a dimcnsionless perturbative parameter which becomes 
relevant (in the RG sense) for K > 1/2 and a = Skp 
[25H32]. Interestingly, these RG equations are analogous 
to those describing the commensurate-incommensurate 
transition (CICT) in doped ID Mott-insulating systems 
after the transformation $ < — !• tp and K < — !• 1/ (2K) 
[33TI35] . The crucial term Spx in Eq. ^ plays the role 
of the particle-doping (relative to half-filling case) in the 
CICT, which has the effect of closing the Mott insulating 
gap. Analogously, in our case a finite Sp may close the 
SC gap. The condition a = Skp (i.e. Sp = 0) determines 
the optimal supercurrent J° p through 



Jr. 



3V3 



1 - {S6k F ) 



(4) 



for which the Majorana-carrying topological phase is 
maximally restored. This result is independent of inter- 
actions, and relies on the linearization of e\. ^ around Ep 
(non-linearities may slightly correct the value of J° p ). 

We next analyze the case of non- vanishing Sp and es- 
timate the critical value 8p c for the topological phase 
transition. At very small Spa(£) <C I, the sin function 
in Eq. ([3| is weakly oscillating and the term Spx can be 



dropped, rendering the RG equations similar to the those 
of the (undoped) sine-Gordon model [29 32 . In that case 
and for K > 1/2, we reach the strong-coupling regime 
y {£*) - I at the scale £* = (2 - if" 1 ) -1 In (£ nw /ao) with 
£nw = u/|Ap|, where the SC term ApSin2 , i9 dominates in 
Eq. @. In this regime, the value of $ (x) is pinned to 
the classical minima •§ (x) — {— 7r/4, 3tt/4} of the sin2i9 
potential, reflecting the underlying Z2 symmetry of the 
Majorana chain in the limit L — > 00 [HJ [3TJ US]- As Sp 
increases, the regime Spa (£) > I is eventually reached 
and the sin function becomes strongly oscillating and av- 
erages to zero. At that point the above RG equations 
are no longer valid and the renormalization of y (£) must 
be stopped [37]. The critical value Sp c can be estimated 
from the condition Sp c a (£*) = 1, which implies that 



Sp c 



1 / «o 

&0 \£nw 



K- 



(5) 



This is an important result in our work. In particu- 
lar, the noninteracting case g = (or K = 1) results 
in Sp c oc oc Ap , which has been confirmed by direct 
numerical calculation in the noninteracting model. In the 
case K 1, and for fixed yo = y (t = 0), we observe that 
repulsive (attractive) e-e interaction destabilizes (stabi- 
lizes) the topological phase, inducing a smaller (larger) 
Sp c . Importantly, for K > 1/2 and tilt-angle \9\ < 7r/2, 
Eq. ([5]) implies that the topological phase can always be 
restored with a supercurrent such that \a — Skp\ < Sp c . 

We consider now the experimental consequences of our 
findings. As mentioned before, the existence of MZBSs 
can be experimentally verified with tunneling transport 
spectroscopy [T7H20] . We consider a single normal metal- 
lic lead with a bias- voltage eVf,, weakly coupled to the left 
end of the NW via the tunneling Hamiltonian 



Hp — 



p>i 



where T is the tunneling coefficient matrix, d p is the elec- 
tron annihilation operator in metallic lead, 7j are the 
second-quantization MZBS operators localized at the left 
(j = L) and right (j = R) ends of the NW. The sum J2' 
runs over the ID-bulk states in the NW, and the coupling 
coefficients \T p x\ ^S> \T p ^r\ due to the exponentially local- 
ized Majorana wave functions. In the topological phase, 
both the MZBS and the ID-bulk continuum modes in 
the NW contribute to the tunnel current /. Using the 
Keldysh formalism, we obtain the tunnel current from 
I = -eN = -f[H T ,N], where N = J2 P d l d P is the 
number of electrons in the metallic lead. Following Refs. 
[TH1 |2"U] we obtain the expression 



I = 



■ J dujTr [r e g R (w) v h g A (w)] [1 - / (w - eV h )\ 



+ - / dujT H N H [1 - / (w - eV b )} , (7) 
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where / (w) is Fermi distribution function and the 
trace is taken in the subspace spanned by jj modes. 
The retarded and advanced Majorana Green's func- 
tions G r {uj) = [£ A (w)] f and [g^H] -1 = oj/2 + 

i [r e (w) + r h ( w )] /2, where r?. ( w ) = rj 1 / (-w) - 

2 7r Sp^p,i^pj^( w — £ p) are the self- energies, and e p the 
single-electron dispersion relation in the metallic lead. 
The second term in the right hand side of Eq. ([T]) 
represents the contribution from ID-bulk states, where 
r(w) = 2irJ2 p \Tp, q \ 2 6{uj - e p ), and N (u) is the ID - 
bulk density of states in the NW. 
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FIG. 2: (Color online) dl/dV for (a) = and (b) 9 = 0.2vr 
with J s = 0. (c,d) Restoring the ZBP at $ = 0.2tt by su- 
percurrents. The blue, red, black, and green curves corre- 
spond to the temperature T = 0, 60mK, 180mK, and 360mK, 
respectively. Other parameters are taken as Vo = l.OmeV, 
E so = O.lmeV, A s = 0.5meV, and the tunneling energies 
|r e r f| ~ irl = 0.005meV. 



Numerical results of dl/dV are plotted in Fig. [2] (a-d) 
at different temperatures. For J s = 0, a ZBP is obtained 
when <6 C ~ 0.1687T (Fig. [2] (a)), and disappears when 
9 > 9 C (Fig. [2] (b)). This result is consistent with the 
observation in the recent experiment by Kouwenhoven 
group [3T]. Fig.[2](c,d) shows that the Majorana-carrying 
phase is restored by a finite supercurrent along +x direc- 
tion at 9 = 0.27T, and maximizes at J s = J° p w 0.039J C 
with £ - lOnm (Fig.|](d)) 25,. The ZBP in the dl/dV 
spectra is clearly restored, indicating the reemergence of 
MZBS after the bulk SC gap reopens. Further increasing 
J s again reduces the bulk gap (see Fig.[I](b)), and finally 
closes it. We confirm that the ZBP is 2e 2 /h at T = 0, 
when the tunneling coefficients are small relative to the 
superconducting bulk gap, and is strongly suppressed by 
thermal broadening. The disappearance and restoration 
of the ZBP provide useful experimental tests for topolog- 
ical superconductivity in the lab. 

Finally we propose an important application of our 
findings to the braiding of MZBS, as needed in TQC. 





FIG. 3: (Color online) Supercurrent assisted braiding of two 
MZBS in a "y" junction, (a) For 9 > 9 C , the NW segment L 2 
is initially in the trivial phase at J s = 0. (b) Move the MZBS 
71 to NW L2 by applying a supercurrent J a — J° p in L 2 . (c) 
Move 72 to the original position of 71. (d) Move 71 to the 
NW Li, and then turn off the supercurrent. 



For a ID system, the braiding operation of MZBS in a 
single NW is not well defined, and the minimum require- 
ment to exchange two MZBS is to consider a "T" or "Y" 
junction composed of several NW segments [351 [37] ■ A re- 
alistic 2D/3D network of MZBS applicable for TQC can 
in principle be constructed by putting together multiple 
NW junctions [38] . However, in such a network there 
must be NW segments which are unavoidably misaligned 
with respect to the external magnetic field, and therefore 
the SIS is generically broken in those NWs. Thus, being 
able to drive all NWs deep into the topological phase 
then becomes questionable, bringing an inevitable diffi- 
culty to braid MZBS in a realistic network. To resolve 
this issue, we introduce a supercurrent-assisted braiding 
(SAB) scheme, as shown in Fig. [3] (a-d) with a "y"- 
junction. In this setup, the quantization axis of eigen- 
spins of a Rashba-type SO coupled NW is perpendicular 
to the NW and parallel to the SC plane (interface of the 
SC/NW heterostructure) [5T]. To minimize the orbital 
effects, the external magnetic field B must lie in the SC 
plane, therefore breaking SIS for at least one of the two 
NW segments. If the B field is applied along the NW seg- 
ment L\ (Fig. |3| , the segment L 2 is topologically trivial 
at J s — when 9 > 9 C . On the other hand, to avoid the 
existence of MZBS and low energy Andreev bound states 
at the junction point of Li,2, the tilt angle 9 should be 
as close to ir/2 as possible [3S]. For the same parameter 
as taken in Fig. [T] for InSb NW [21], we have the critical 
angle 9 C w 0.168tt. Then for = 0.2tt, the NW L 2 is 
already in the trivial phase without applying a supercur- 
rent (Fig. [3]ja)), and we next exchange two MZBS 71,2 
localized on the ends of L\. To perform the braiding of 
71.2, we turn on a supercurrent with J s = J° p ps 0.039J c 
along L2 (for £ ~ lOnm in NbTi s-wave SCs [25]) and 
move first 7! to NW L 2 by gate control (Fig. [3] (b)). 
Then we move 72 to the original position of 71 (Fig. [3] 
(c)). Finally 71 is shuttled to L\, completing exchange 
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with 72, and the supercurrent is turned off after braid- 
ing (Fig. (3jd)). It is noteworthy that the supercurrent is 
needed only in the intermediate process of the braiding 
operation. Applying the SAB to generic 2D or 3D Majo- 
rana networks can provide vast flexibility for the realistic 
TQC with MZBS. 

In summary, we have studied the disappearance and 
reemergence of MZBS in Majorana quantum wires with 
broken SIS due to the application of a tilted external 
field and supercurrents. We have shown the robustness 
of these findings against the presence of e-e interactions, 
where the physics is analogous to the case of doped ID 
Mott-insulating systems. The present study provides 
new insights into the braiding operation in a realistic net- 
work of MZBS, where the SIS is generically broken. We 
introduced the supercurrent assisted braiding of MZBS, 
and show that it has important applications to the real- 
istic Majorana- fermion-based quantum computation. 

We thank Meng Cheng for helpful comments. We 
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